Red Density Perturbations and Inflationary Gravitational Waves 
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We study the implications of recent indications for a red spectrum of primordial density per- 
turbations for the detection of inflationary gravitational waves (IGWs) with forthcoming cosmic 
microwave background experiments. We find that if inflation occurs with a single field with an 
inflaton potential minimized at V = 0, then Planck will be able to detect IGWs at better than 2er 
confidence level, unless the inflaton potential is a power law with a very weak power. The proposed 
satellite missions of the Cosmic Vision and Inflation Probe programs will be able to detect IGWs 
from all the models we have surveyed at better than 5<r confidence level. We provide an example of 
what is required if the IGW background is to remain undetected even by these latter experiments. 



Inflation's [l[ predictions of a flat Universe and a 
nearly Gaussian and nearly scale-invariant spectrum 
of primordial density perturbations [2j have been con- 
firmed by a suite of recent cosmic microwave back- 
ground (CMB) experiments @, 0, [1]. A third pre- 
diction of inflation is a stochastic background of in- 
flationary gravitational waves (IGWs) |7j- The the- 
ory predicts that the IGW amplitude is proportional 
to the cosmological energy density during inflation. 
Since this energy density varies from model to model, 
the IGW amplitude cannot be predicted in a model- 
independent way. If, however, IGWs can be detected, 
they will provide an important probe of the physics of 
inflation. 

The CMB polarization is perhaps the most promis- 
ing tool for detecting the IGW background. The sta- 
tistical properties of CMB linear polarization are fully 
characterized by two sets of spin-2 multipole moments 
with opposite parities Q. The magnetic-type modes 
(B or curl modes) are produced by IGWs (or "tensor" 
metric perturbations) and not density perturbations 
(scalar metric perturbations), and they do not corre- 
late with the temperature nor the electric-type-parity 
modes (E or grad modes) . A detection of B-mode po- 
larization would thus provide good evidence for IGWs 
9]. 

The scalar and tensor perturbations generated by 
inflation have power spectra that are generally well ap- 
proximated by power laws: P s (k) oc k n " and P t (k) oc 
k nt , respectively, as a function of the spatial wavenum- 
ber k. The amplitudes and spectral indices can be 
written in terms of the inflaton potential V(cj>) and its 
first and second derivatives V'(4>) and V" ((/))■ More 
precisely, the spectral indices n s and n t can be written 
in terms of the slow-roll parameters e and r\ 10], 
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where mpi is the reduced Planck mass. The observed 
density-perturbation amplitude already fixes V/e ~ 
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6.6 x 10 16 GeV. The tensor-to-scalar ratio r is r = 
16e. Since n% and r depend only on e, they satisfy 
a consistency relation n t ~ — 2e = — r/8 that can be 
tested with measurements, if the IGW background can 
be detected. 

An important new result from recent CMB exper- 
iments is that a value n s — 1 now seems unlikely; 
the likelihood for its value peaks closer to n s — 0.95 
UGH. Although it is still debated whether the data 
are inconsistent with n s = 1 at the 2a level or at the 
3(7 level [12|, it is clear that a value n s = 1 is becom- 
ing increasingly difficult to reconcile with the data. 
Although a measurement of n s does not uniquely de- 
termine e, it is natural to expect that e and 77 are 
roughly of the same magnitude, or perhaps that e 
is much larger than 77 (if, e.g., the inflaton potential 
can be approximated as linear). Conversely, it would 
be unusual to have a value rj 3> e, as it would re- 
quire that the inflaton <j) sit near an inflection point 
of V(4>), an unusual location. A value of <f> near an 
inflection point is also unworkable in many potentials, 
as it does not allow the right number N e of e-foldings 
between CMB scales and the end of inflation (which 
occurs when the potential becomes sufficiently steep). 
If e > rj (and barring some unusual cancellation) , then 
a value n s ~ 0.95 - 0.99 implies e = 0(0.01), which 
thus implies a r ~ 0.1, perhaps within the range of 
detectability with next-generation experiments, as ar- 
gued similarly in Ref. [13| . 

The purpose of this paper is to quantify these ar- 
guments a bit more precisely by considering a vari- 
ety of functional forms for V((f>), consistent with cur- 
rent data. As we will see, the amplitude of the IGW 
background is large enough to be detectable by next- 
generation experiments in a broad family of inflation 
models consistent with current data. 

We consider three classes of inflationary models 
0, [Ulillij]: (1) Power-law inflation is characterized 
by an inflaton potential V{4>) oc e^/ M , where /1 is 
a mass scale. In this potential, there is a relation, 
r = 8(1 — n s ), between the tensor-to-scalar ratio and 
the scalar spectral index. (2) Chaotic inflation fea- 
tures an inflaton potential V(<f>) oc (0//i) p . In theo- 
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retically attractive models, p is a small integer. Ex- 
perimentally, p < 10 if n s > 0.9 [l6j], but p can em- 
pirically be arbitrarily small. We consider values of p 
between p = 0.1 and p = 8 in our numerical work. In 
these models, r — 8[p/(p + 2)](1 — n s ). (3) Sponta- 
neous symmetry-breaking (SSB) inflation features an 
inflaton potential V(4>) oc [1 — (4>/v) 2 ] 2 . The precise 
model is specified by two parameters: v and iV e , the 
number of e-foldings of inflation between CMB scales 
and the end of inflation. A conservative range for N e 
is 47 < N e < 62, corresponding to an inflationary en- 
ergy scale in the range of an MeV to 10 16 GeV, the 
current upper bound. The n s -r relation for SSB mod- 
els cannot be written in a simple way. To obtain it, we 
use the algorithm given in Eqs. (38)-(41) of Ref. jig . 
We show in Fig. [1] the corresponding regions in the 
n s -r parameter space. 

One could also consider hybrid-inflation models, 
which generally feature multiple fields. Their phe- 
nomenology can usually be modeled, though, as a sin- 
gle field with the addition of a non-zero cosmological 
constant at the inflaton-potential minimum. What 
distinguishes these models phenomenologically is that 
they usually produce n s > 1. Since the working as- 
sumption of this paper is that n s < 1, we do not 
consider these models further, but simply note that 
if n s is indeed greater than unity, than these mod- 
els allow for much smaller IGW amplitudes than the 
single-field models we consider here. 

In their analysis Q of their three-year data, the 
Wilkinson Microwave Anisotropy Probe (WMAP) col- 
laboration surveyed the cosmological parameter space 
assuming n s , r, and nt to be independent parameters. 
The resulting constraints to the n s -r parameter space 
are shown in Fig. [T] However, if inflation is respon- 
sible for primordial perturbations, and if the inflaton 
potential is power-law, chaotic, or SSB, then there is 
a fixed relation between n s , r, and n*; i.e, they are 
not independent parameters. 

We therefore re-analyze the CMB data for the 
power-law, chaotic, and SSB models, imposing the 
consistency relations between n s , r, and nt- To do 
so, we use the Markov chain Monte Carlo (MCMC) 



package cosmomcmc 17J to run a set of chains, impos- 
ing these consistency relations. We use only the lat- 
est WMAP results [4(. The likelihood is determined 
using the October 2006 version of the WMAP like- 
lihood code available at Ref. [5|. The likelihood are 
obtained after marginalizing (with flat priors) over the 
baryon and cold-dark- matter densities, the ratio of the 
sound horizon to the angular-diameter distance at de- 
coupling, and the optical depth to reionization. 

In TablelH we list the current constraints to n s and r 
from WMAP obtained using each class of inflationary 
models as a prior. As expected, the tabulated values 
and the plots indicate that tensor modes are required 
for exponential, power-law, and small-field scenarios. 

Our analysis provides a likelihood, from CMB data, 
for each point along the curves associated with each 




FIG. 1: The WMAP (red shaded areas) constraints to 
the n s -r parameter space. The outer (inner) contours are 
the 2a (la) contours. We also plot the models within 
A\ 2 = 1 (lo") from the best fit for power-law and chaotic 
models, and at A\ 2 = 2.3 (also la) for the SSB models. 
The curves, from top to bottom, are for power-law and for 
chaotic p = 8, p = 1, and p — 0.1. The (green) points are 
for the SSB models. 
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FIG. 2: Percentage of inflationary models with a tensor- 
to-scalar ratio above a threshold value r, for the various 
inflaton-potential classes considered in Fig. [1] From right 
to left, the models are power law, chaotic p = 8, chaotic 
p = 1, SSB, and chaotic p = 0.1. 



class of inflationary models. We than thus plot in 
Fig. [2] the percentage of models with a tensor ampli- 
tude above a threshold value r. As shown, almost all 
of the exponential and power-law models have a ten- 
sor amplitude r > 2 x 10~ 2 , while small-field models 
predict r > 2 x 10" 3 . If experiments were to probe 
values of r well below 2 x 10~ 3 without finding any 
evidence for gravitational waves, then it would rule 
out a large class of single-field inflationary models. 
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Model 


n S 


r 


Power Law 


0.980_ 005 




Chaotic p=l 


0.970±g;SS| 


0-080lg;SS 


Chaotic p=8 




14+ ' 07 


Chaotic p=0.1 


964+ - 019 


014+ ' 004 
u,ul *-0 .007 


SSB (A^ e = 47 - 62) 


U M0( -0.020 


u.u4z_ 014 



TABLE I: The 68% C.L. limits on the spectral index n s 
and the scalar-tensor ratio r from WMAP, assuming the 
n a -r relation for each class of inflationary models. 



Previous work [18J can already be used to estimate 
the fraction of models in each class of models that can 
be detected by a given experiment. Here we go further 
by considering values of rit ^ that are given by the 
inflationary consistency relation. We also evaluate the 
likelihood, given current WMAP constraints to each 
class of inflationary models (e.g., power-law, chaotic, 
SSB), that IGWs will be detected by a particular fu- 
ture experiment. 

We begin by evaluating, for a given future experi- 
ment and for a given set of inflationary parameters, 
the corresponding Fisher Matrix [l9( 

«^EEf«;f. (3) 

1=2 a, 13 1 3 

where the C" 13 are the power spectra for the temper- 
ature (TT), temperature-polarization (TE), E-mode 
polarization (EE), and B-mode polarization (a and 
run over TT, EE, TE, BB), 6i is the set of parameters 
used in the Markov Chain, and Cov/ is the "spectra 
covariance matrix" 0] . The covariance-matrix entries 
depend on the instrumental noise and the angular res- 
olution of the experiment. The la (or 68% C.L.) error 

on a particular parameter 6j is cr(6j) — {F^ 1 ) 1 /^ [Hj], 
after marginalizing over other undetermined parame- 
ters. 

Now consider a particular class (e.g., power-law, 
chaotic, or SSB) of inflationary models. The prob- 
ability that IGWS will be detected at Na, given the 
current WMAP constraints, is then, 

M 

P = Y,®[n-N<j{n)}/M, (4) 

i=l 

where © is the Heaviside step function. The sum 
is over the M models consistent at 95% C.L. with 
WMAP assuming one of class of inflationary models. 
For example, if P = 1 for N = 2, then the future ex- 
periment will detect IGWs with 2a significance from 
all models compatible with WMAP within that par- 
ticular class of models. 

We consider here four experimental configurations: 
eight years of WMAP; the Planck satellite [20] (to be 
launched in 2008); and two proposed next-generation 



Experiment Chan. FWHM AT/T AP/T 



WMAP 


40 


28' 


8.2 


11.6 


/ sky = 0.65 


60 


21' 


11.0 


15.6 




90 


13' 


18.3 


25.9 


Planck 


44 


23' 


2.7 


3.9 


/ sk y = 0.65 


70 


11' 


4.7 


4.7 




100 


9.5' 


2.5 


4.0 




143 


7.1' 


2.2 


4.2 




217 


5.0' 


4.8 


9.8 




353 


5.0' 


14.7 


29.8 


BPol 


45 


900' 


0.013 


0.026 


/sky = 0.65 


70 


68' 


0.12 


0.24 




100 


47' 


0.062 


0.12 




143 


47' 


0.055 


0.11 




217 


40' 


0.046 


0.092 




353 


59' 


0.062 


0.12 


EPIC 


40 


116' 


0.032 


0.047 


/sky = 0.65 


60 


77' 


0.018 


0.039 




90 


52' 


0.0077 


0.025 




135 


34' 


0.0073 


0.036 



TABLE II: Experimental specifications for the satellite 
missions considered in this work. Channel frequency is 
given in GHz, FWHM in arcminutes, and noise in 10~ 6 . 



satellites, BPol [2l| and EPIC [22j. Since the ampli- 
tude of the Galactic foregrounds for polarization is 
still not determined by observations, we remove chan- 
nels below 40 GHz and above 250 GHz, as these are 
likely to be contaminated by synchrotron and dust 
emission, respectively. The parameters assumed for 
these missions are summarized in Table [TT1 



Model 


WMAP8yr 


Planck 6ch BPol Epic 


Power-law 


0.45/0 


1/0.98 1/1 1/1 


Chaotic p=l 


0/0 


0.99/0.67 1/1 1/1 


Chaotic p=8 


0.30/0 


1/0.97 1/1 1/1 


Chaotic p=0.1 


0/0 


0.60/0 1/1 1/1 


SSB (N e = 47 - 62) 


0/0 


0.78/0.09 1/1 1/1 



TABLE III: Percentage of models in agreement with the 
WMAP observations and with an IGW background de- 
tectable at 2<j and 5a confidence levels by the experimental 
configurations listed in Table ITT1 



Our results are reported in Table IIIII for detec- 
tions with 2a and 5<r significance. Our results sug- 
gest that eight years of WMAP data are unlikely 
to detect IGWs. The Planck satellite, on the other 
hand, stands a good chance to measure the IGW 
background in power-law models, even at the 5a con- 
fidence level. This forecast, however, may be opti- 
mistic, since it assumes perfect foreground removal 
and full use of 6 channels for cosmology. To under- 
stand these effects better, wc show forecasts in Ta- 
ble IIVI for Planck, assuming first a smaller number 
of channels and then a more pessimistic configuration 
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iviocei 


Ar-h 


lUUvjXlZCIl 


alien _l -\- 


Power Law 


1/0.97 


0.99/0.30 


0.81/0.02 


Chaotic p=l 


0.99/0.61 


0.92/0 


0.25/0 


Chaotic p=8 


1/0.96 


0.99/0.16 


0.81/0 


Chaotic p=0.1 


0.50/0 


0/0 


0/0 


SSB (7V e = 62 - 47) 


0.77/0.07 


0.36/0 


0/0 



TABLE IV: Planck Only. Percentage of models in agree- 
ment with the WMAP observations and with a GW- 
background detectable at 2<r/5 a by the Planck experiment 
in function of the channels available for cosmology. The 
T + C does not include the polarization channels but just 
the temperature and the cross temperature-polarization 
spectrum. 



just using TT and TE spectra alone, will not be able 
to produce more than an indication for tensor modes. 
With either an Inflation Probe or a Cosmic Vision 
mission for high-precision CMB polarization measure- 
ments, we find that the full set of single-field inflation- 
ary models we have considered can be explored. If no 
tensor modes are detected after one of these missions, 
then it is very likely that inflation cannot be described 
by a single field that rolls to a minimum with zero 
cosmological constant. The alternative would then 
be more complicated models, such as chaotic inflation 
with an unusually small p (p < 0.01), models with 
multiple fields, string models, or models with extra 
dimensions. 



where no EE or BB power spectra are used. While we 
found only a small degradation in decreasing the num- 
ber of channels from 6 to 4, using one channel only or 
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